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Abstract 


■-  v  n  n 

— »  Kinetic  modeling  calculations  on  the  H2-02  system  have  been  carried  out 
with  an  extensive  reaction  set  to  probe  the  vicinity  of  the  three  explosion 
limits.  Sensitivity  analysis  is  used  throughout  this  investigation  to 
study  system  behavior,  in  particular,  to  elucidate  mechanistic  details.  The 
concentrations  and  sensitivity  profiles  are  discussed  in  light  of  the 
appropriate  experimental  results  and  existing  theories  of  hydrogen  combustion. 
The  results  indicate  the  present  model  to  be  useful  over  a  wide  pressure- 
temperature  range i^-The  reaction  set  is  also  used  to  probe  the  sensitivities 
for  an  experimental  study  designed  to  measure  the  rate  constant  of  an 
important  elementary  reaction,  H  +  02  +  M  — >  H02  +  M  ,  involved  in  this 
system.  The  versatility  of  the  reaction  set  is  also  demonstrated  by  a  study 
of  a  related  chemical  reaction,  the  decomposition  of  hydrogen  peroxide. 
Finally,  prospects  for  utilizing  the  methods  and  results  of  this  study  to 
examine  other  complex  kinetic  schemes  are  discussed. 


I.  Introduction 


Numerical  modeling  is  becoming  a  valuable  tool  in  elucidating  the 
behavior  of  large  and  complex  systems  in  fields  as  diverse  as  reaction 
kinetics,1  operations  research,2  world  models  for  population  and  resources,3 
and  mathematical  biology.4  Modeling  has  been  shown  to  be  of  particular 
Importance  in  investigating  the  details  of  complex  combustion  processes.5  In 
this  work  we  present  a  model  developed  to  investigate  the  kinetics  of  a 
simple  combustion  process, 

H2(g)  +  *02(g)  - >  H20(g) 

He  present  the  results  of  detailed  numerical  calculations  Involving  the 
proposed  mechanism  for  this  system  over  a  wide  range  of  temperatures  and 
pressures.  These  calculations  contribute  to  an  understanding  of  this  well* 
studied  intriguing  combustion  reaction. 

This  study  of  the  H2-02  system  differs  in  certain  respects  from  similar 
numerical  studies  on  complex  kinetic  schemes6  because  sensitivity  analysis 
plays  a  dominant  role.7  In  any  numerical  study  it  is  desirable  to  understand 
how  sensitive  the  system  solutions  are  to  the  various  input  parameters  of  the 
model.  Generally  most  numerical  studies  of  complex  systems  probe  sensitivities 
via  the  "brute  force"  technique  of  changing  parameters  and  observing  the 
effects  upon  the  output  solutions.  However,  recently  in  the  field  of  chemical 
kinetics,  several  more  systematic  and  economical  approaches  to  sensitivity 
analysis  have  been  developed.8"14  In  particular,  the  Green's  Function  Method 
(GFM)13  has  already  been1 shown  to  be  a  useful  and  competitive  approach  to  the 
problem  of  sensitivity  analysis  in  chemical  kinetics14  as  well  as  in  quantum 


dynamics.15  The  present  paper  provides  an  application  of  sensitivity  analysis 
utilizing  the  GFM  to  deduce  information  regarding  explosion  limits,  mechanistic 
details,  etc.,  for  the  H2-02  system. 

In  brief,  then,  this  work  investigates  what  we  believe  to  be  a  realistic  * 
model  for  hydrogen  combustion  by  means  of  a  detailed  analysis  of  concentrations 
and  sensitivity  profiles.  An  attempt  is  made  to  identify  and  understand  the 
somewhat  enigmatic  temperature-pressure  behavior  of  this  system,  the  Important 
mechanistic  pathways  and  the  critical  rate  constants  for  the  process.  The 
numerical  results  -  concentrations  and  sensitivities  -  are  discussed  in  light 
of  relevant  experiments  and  models.  The  utility  of  sensitivity  analysis  and 
kinetic  modeling  is  apparent  by  suggestions  from  this  work  for  additional 
experiments  to  understand  this  important  system  more  thoroughly.  Finally, 
the  reaction  set  developed  for  the  H2-02  system  is  shown  to  be  versatile 
enough  to  analyze  sensitivities  in  an  experimental  study  of  an  Important 
elementary  step  H  +  02  +  M  — *  H02  +  M  of  the  model.  Calculations  have 
also  been  carried  out  to^ probe  a  related  chemical  reaction,  the  decomposition 
of  hydrogen  peroxide. 

The  paper  is  divided  into  six  remaining  sections.  Section  II  provides  a 
brief  introduction  to  the  H2-02  system  and  discusses  the  development  of  the 
reaction  model.  Section  III  discusses  some  computational  details  for  this 
system.  Section  IV  is  a  detailed  analysis  of  the  results  of  many  calculations, 
focusing  upon  system  behavior  at  the  three  explosion  limits  for  the  reaction. 
Section  V  is  an  analysis  of  sensitivities  for  an  experimental  probe  of  the 
H  +  02  +  H  - *  H02  +  M  step  in  the  model.  Section  VI  discusses  the  results 
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of  calculations  for  the  decomposition  of  H202.  Finally,  Section  VII  provides 


II.  An  Introduction  to  the  Ho-Oo  System 


The  formation  of  HgO  from  Hg  and  02  has  been  the  subject  of  hundreds  of 
papers  and  is  still  a  problem  for  active  research.  At  room  temperatures  and 
ordinary  pressures,  the  reaction  is  very  slow.  When  the  mixture  is  heated 
and  subjected  to  a  spark  or  an  activating  surface,  however,  the  process  can 
be  very  rapid,  even  explosive.  Particularly  interesting  behavior  occurs  in 
the  temperature  range  650  -  850°K  for  reactant  pressures  of  1  to  5,000  torr. 
Figures  1  and  2  illustrate  some  peculiarities  in  the  behavior  of  a 
stoichiometric  mixture  of  H 2  and  Og.  The  so-called  first,  second  and  third 
explosion  limits  are  shown.  Explosion  limits  are  not,  strictly  speaking, 
flamnability  limits.  They  are  the  pressure-temperature  boundaries  for  a 
specific  mixture  of  fuel  and  oxidizer  that  separate  the  regions  of  slow  and 
fast  reactions.13  In  the  explosive  regime  flames  can  be  supported  via 
steady  diffusional  mechanisms.  Many  experiments  have  been  carried  out  to 
probe  the  effects  of  vessel  size,  vessel  composition,  surface  composition, 
inert  gases,  etc.,  on  these  three  limits.  Although  these  experiments  are 
described  in  considerable  detail  in  several  texts, 17-22  a  brief  summary  of 
the  results  for  each  explosion  limit  is  appropriate. 

Below  the  first  explosion  limit  the  reaction  proceeds  almost  negligibly 
slowly.  However,  if  a  stoichiometric  mixture  of  Hg  and  02  is  added  slowly  to 
a  reaction  vessel,  the  mixture  experiences  a  mild  explosion  near  a  few 
torr.  The  pressure  is  defined  as  the  sum  of  the  initial  partial  pressures 
of  Hg  and  02.  The  first  limit  boundaries  between  the  regions  of 
8 low  and  fast  reaction,  is  quite  dependent  upon  vessel  size  and  vessel 
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composition.  Frost  and  Alyea23  found  that  rinsing  their  vessels  with  KC1 


raised  the  limit  pressure  to  several  times  higher  than  previously.  The 
marked  dependence  upon  vessel  size  and  vessel  composition  indicates  that 
diffusional  mechanisms  are  very  important  in  this  regime.  There  is 
competition  between  the  formation  of  radicals  which  undergo  chain  reactions 
and  the  destruction  of  such  chains  via  migration  to  the  walls  of  the  vessel. 

The  neighborhood  of  the  second  explosion  limit  has  been  probed  by  both 
the  heating  method,  in  which  a  mixture  of  H2  and  02,  initially  in  the 
nonexplosive  regime  between  the  second  and  third  limits,  is  heated  until  an 
explosion  is  observed,  and  by  the  withdrawal  method,  in  which  slow  evacuation 
of  the  mixture  is  carried  out  until  explosion  occurs.21  The  experimental 
evidence  indicates  that  behavior  in  the  neighborhood  of  the  second  explosion 
limit  must  be  explained  by  gas  phase  production  and  destruction  of  radicals.11 
Evidently,  diffusion  plays  a  less  dominant  role  than  it  does  for  the  first 
explosion  limit,  since  the  second  limit  pressure  is  independent  of  the 
vessel,  size,  provided  its  volume  is  sufficiently  large.17  Still,  there 
is  some  dependence  upon  the  composition  of  the  surface,  that  is,  different 
kinetic  behavior  is  observed  for  uncoated  surfaces  and  those  coated  with 
salt  or  acid.17  This  indicates  the  important  role  played  by  the  surface 
in  the  initiation  and,  to  a  lesser  extent,  destruction  of  gas  phase  radicals. 
Once  produced,  the  production  of  these  radicals  is  accelerated  by  chain 
branching  reactions  such  as  H  +  02  — »  OH  +  0  ,  resulting  in  very  rapid 
reaction.  The  slow  reaction  at  higher  pressures  is  usually  explained  by 
a  third  order  non-branching  reaction  which  overtakes  the  H  +  02  — *  OH  +  0 
step.16  Although  some  Investigators  have  attempted  to  describe  the  second 
limit  without  recourse  to  the  metastable  H02  radical21/24,  it  is  now 


generally  accepted  that  the  important  third  order  reaction  is 

H  +  02  +  M  - >  H02  +  M 

The  characteristic  slow  reaction  at  pressures  between  the  second  and  third 
limits  proceeds  through  the  formation  of  the  hydroperoxyl  radical  and 
hydrogen  peroxide,  which  has  been  observed  experimentally  in  this  regime. 25 

At  the  first  and  second  limits,  the  observed  explosions  are  of  the 

isothermal  branched-chain  type,  and  just  outside  the  limits  the  overall 

reaction  rate  is  extremely  small.21  At  the  third  limit,  however,  there  is 

no  such  sharp  transition  from  slow  reaction  to  explosion,  and  the  reaction 

rate  is  high  even  below  the  limit.17  For  this  reason  it  has  not  been 

established  whether  the  third  limit  has  mainly  a  chain  or  a  thermal  character 

Potassium  chloride  wall  coating  hinders  the  explosion,  suggesting  that  chains 

are  involved.  Careful  studies  of  the  third  limit  are  difficult  because  of 

technical  problems  in  handling  the  reactants  and  determining  their  pressures 

prior  to  the  onset  of  a  very  rapid  reaction.  The  explosion  is  thought  to  be 

due  to  a  combination  of  chain  and  thermal  mechanisms,  that  is  to  say,  the 

heat  given  out  by  the  chain  reactions  cannot  be  dissipated  sufficiently 

rapidly,  so  that  self-heating  occurs,  followed  by  explosion.22  Glassman16 

suggests  that  the  third  limit  is  due  to  a  reaction  that  overtakes  the 

stability  of  the  H02,  most  possibly  the  reaction 
* 

HOg  +  Hg  HgOg  +  H 

which  regenerates  a  large  number  of  active  particles.  Lewis  and  Von  Elbe21 
have  carried  out  some  experiments  near  the  curve  joining  the  second  and  third 


explosion  limits.  One  experiment  investigates  a  stoichiometric  mixture  of 
H2  and  02  in  N2  initially  at  847°K  and  J>00  torr,  for  which  the  initial 
reaction  rate  (rate  when  ~  10$  of  the  H20  has  been  formed)  is  about 
1.5  torr/min.  In  this  case  no  explosion  occurs,  as  the  rate  is  too  small  to  - 
seriously  disturb  isothermal  conditions:  the  heat  is  rapidly  dissipated. 
However,  at  823° K  and  950  torr,  the  initial  rate  is  20  -  40  torr/min,  which 
is  sufficiently  large  to  raise  the  temperature  and  then  cause  explosion 
after  a  very  short  induction  period.  This  particular  experiment  is  described 
in  considerable  detail,  since  very  reliable  experimental  results  are  lacking 
in  this  regime  because  of  the  difficulties  we  have  cited.  The  computational 
results  described  in  Section  IV  will  shed  some  light  on  kinetic  behavior  in 
this  regime. 

There  have  been  many  attempts  to  explain  the  complicated  pressure- 
temperature  behavior  of  this  system  by  various  models .  Some  models  include 
very  few  reactions  and  attempt  to  include  diffusional  mechanisms,  heat 
release,  etc.28  Others  have  focused  primarily  upon  mass  action  kinetics. 
Baldwin,  et  al.,27  have  constructed  a  reaction  set  which  is  representative 
of  the  types  of  kinetic  mechanisms  that  have  been  considered. 

In  building  the  reaction  set  for  this  paper,  we  have  assumed  a  conserva¬ 
tive  posture  by  including  virtually  every  conceivable  reaction  which  might  be 
important  that  involves  the  species  H2,  02,  H20,  H202,  H02,  H,  0,  and  OH.  All 
these  species  have  been  positively  identified  in  hydrogen  flames,  although 
reactions  involving  ozone  are  often  excluded,  since  its  concentration  rarely 
exceeds  a  few  ppm.26  The  present  authors  have  carried  out  an  extensive 
literature  search  of  several  data  bases  to  construct  the  reaction  set  and  to 


obtain  the  best  possible  numerical  values  available  for  the  rate  constant 
parameters  for  the  elementary  steps  of  the  mechanism.  The  literature  values 
have  been  carefully  reviewed  and  values  assigned  based  upon  the  authors' 
evaluations  regarding  which  measurements  and/or  estimates  are  the  most 
accurate.  The  results  of  this  effort  are  embodied  in  Table  I,  which  lists 
each  reaction  with  appropriate  rate  constant  parameters  and  uncertainties. 

The  reaction  set  of  Baldwin,  et  al..27  is  also  indicated  in  this  table.  In 
general  the  rate  constants  are  of  a  modified  Arrhenius  form,  i.e., 

Bi 

k£  =  AiT  exp(-Ct/T)  (2.1) 

Most  of  the  calculations  described  in  the  following  sections  of  this  paper 
utilize  the  first  57  reactions  of  Table  I.  The  wall  reactions  (58-62)  have 
been  included  in  the  mechanism  to  probe  wall  effects  in  the  vicinity  of  the 
first  and  second  limits.  Of  course,  the  numerical  values  of  these  wall 
reaction  rate  constants  varies  with  the  vessel  surface- to-volume  ratio,  its 
composition,  the  method  of  preparation,  etc.  Thus  we  have  chosen  values  from 
the  literature64* 85  corresponding  to  values  for  typical  laboratory  conditions 
to  merely  simulate  wall  effects. 

Several  adiabatic  calculations  were  carried  out  in  the  region  of  the 

third  (thermal)  explosion  limit.  These  calculations  required  the  specific 

enthalpy  (h^)  and  the  specific  heat  capacity  (Cv  )  for  each  species.  For  the 

purposes  of  these  calculations  both  h.  and  C  are  adequately  represented 
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as  linear  functions  of  temperature  in  the  range  of  interest  ( 500  -,6000°K) . 
Linear  least  squares  fits  to  the  tabulated  JANAF  data66  were  used  to  obtain 
the  coefficients  listed  in  Table  II. 


The  purpose  of  this  section  is  to  explain  the  nurseries  of  a  single 
calculation  in  considerable  detail  in  order  to  show  how  the  computations 


are  carried  out,  what  sort  of  information  is  obtained,  what  error  tolerances 
are  used,  etc.  Detailed  discussions  of  kinetic  and  sensitivity  results  are 
given  in  Sections  I V,  V  and  VI. 

The  formulation  and  numerical  considerations  of  the  Green's  Function 


Method  of  sensitivity  analysis  have  already  been  treated  in  two  other 


papers. 13 ’ 14  Briefly,  to  determine  concentrations  and  sensitivity  profiles, 
one  must  first  solve  the  initial  value  problem  resulting  from  the  application 
of  the  law  of  mass  action67  to  reaction  mechanisms.  The  resulting  (generally 
stiff)  differential  equations  are  of  the  form 


•hi  =  £(:&  9b  t)  >  z(o)  »  r 


(5.1) 


where  jr  denotes  an  n-vector  of  species  concentrations,  a  is  an  m-vector  of 

system  parameters  (i.e.,  rate  constants  and  initial  concentrations),  t  refers 

to  time  and  f  is  a  nonlinear  function  of  the  concentrations.  Once  the  n 

concentration  profiles  have  been  obtained,  the  sensitivities  are  determined 

by  solving  n  sets  of  n-dimensional  ordinary  differential  equations  for  the 

adjoint  Green's  function  matrix.  In  the  course  of  solving  these  equations, 

the  mesh  chosen  automatically  by  the  differential  equation  solver68  is  used 

as  a  guide  to  construct  the  grid  required  to  evaluate  the  integrals  for  the 
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sensitivity  coefficients  -r—  (t)  .  The  normalized  sensitivity  coefficients 
d/n  y . 


dinar 


-(t)  are  especially  useful  measures  of  how  sensitive  a  particular 


species  concentration  is  to  a  particular  parameter.  In  subsequent  sections 
we  will  explain  how  these  sensitivities  may  be  used  to  study  system  behavior. 

The  sample  calculation  described  here  has  been  carried  out  at  a  constant 
temperature  of  700°K.  The  initial  pressure  is  320  torr,  the  reactants  H2 
and  02  are  initially  in  a  stoichiometric  ratio  and  no  diluent  is  present. 

This  situation  corresponds  to  the  nonexplosive  regime  between  the  second  and 
third  limits.  Only  the  first  57  reactions  of  Table  I  are  used  in  this 
calculation.  All  calculations  described  in  this  paper  are  performed  in 
double  precision  on  an  IBM  36O/9I  computer.  The  error  tolerance  £  =  10“3 
was  used  in  solving  the  differential  equations,  although  occasionally  a 
stricter  tolerance  was  used  to  check  precision. 

The  differential  equations  for  the  concentrations  were  solved  first,  and 
the  results  are  displayed  graphically  in  Figure  3*  This  behavior  is  fairly 
typical  of  the  qualitative  features  expected  and  observed  for  the  slow  regime 
between  the  second  and  third  limits.  This  regime  will  be  considered  more 
fully  in  Section  IV.  The  differential  equations  for  this  problem  are  very 
stiff;  the  stiffness  ratio  -  the  ratio  of  the  largest  (in  magnitude)  to 
smallest  eigenvalue  -  is  typically  ~101°-1011,  yet  still  only  a  few  hundred 
steps  are  required  by  the  ODE  solver  to  reach  9 0$  reaction  completion.  In 
the  course  of  computing  the  concentrations,  various  terms  contributing  to  f ^  , 
which  are  of  the  form  ±  h^y^y^,  are  computed  and  printed  at  selected 
times.  This  is  useful  in  screening  analysis, 2S  which  assigns  the  importance 
of  reaction  rate  constants  based  upon  the  magnitude  of  the  individual  reaction 
contributions  to  total  species  derivatives.  Mass  conservation  relations14 
are  also  checked  and  verified  at  several  times. 


For  the  sensitivities,  the  differential  equations  for  each  row  of  the 
adjoint  Green's  function  matrix  have  to  be  solved.  This  requires  an 
interpolation  grid  for  the  concentration  vector  £.  Relatively  coarse  grids 
( ~50  grid  points)  are  sufficient  to  obtain  reasonable  precision.  Each  of 
the  n  rows  of  the  adjoint  Green's  function  matrix  -  corresponding  to  the  n 
species  present  -  require  as  many  as  several  thousand  or  as  few  as  twenty 
steps  by  the  differential  equation  solver  to  obtain  sensitivities  at  a 
particular  time.  Sensitivities  are  obtained  at  four  separate  times.  The 
stiff  nature  of  this  system  is  exhibited  in  the  boundary  layer  behavior  for 
the  typical  sensitivity  integrals.  This  type  of  behavior  has  already  been 
documented  for  a  previous  study.14  As  Figure  4  shows,  the  Integral  is 
generally  determined  by  only  a  very  small  portion  of  the  integrand.  Since 
the  exponential  option  for  numerical  quadrature14  is  employed,  and  since  the 
mesh  chosen  by  the  solver  is  used  to  construct  the  Integration  grid  (in  this 
study  every  third  mesh  point  is  utilized),  the  results  should  be  quite 
accurate.  The  sensitivity  mass  conservation  relations  described  in 
reference  14  are  reasonably  well  satisfied. 

The  final  point  concerns  computational  time  estimates.  To  obtain 
concentrations  at  a  reasonably  fine  grid  of  times  for  up  to  90$  reaction 
completion  requires  about  three  seconds  of  execution  time  on  the  IBM  360/91 
computer.  The  adiabatic  calculations  described  later  in  the  paper  require 
about  twenty  seconds.  To  obtain  sensitivities  for  a  single  species  with 
respect  to  all  parameters  at  a  single  time  requires  only  a  few  more  seconds, 
depending  upon  which  species.  The  times  are  essentially  independent  of  the 
number  of  parameters  considered.  Generally,  it  was  found  desirable  to  obtain 


sensitivities  for  all  species  at  least  four  tines  during  the  reaction  for 
this  and  most  of  the  other  calculations.  Since  occasionally  the  solver  takes 
thousands  of  steps  for  one  species,  the  total  e  ?cution  tine  is  about  2.1 
minutes,  which  is  fairly  typical  for  most  of  the  calculations.  Based  upon 
previous  studies,14  it  is  expected  that  the  direct  method  of  sensitivity 
analysis8"10  would  take  at  least  twenty  minutes  of  execution  time  to  obtain 
similar  sensitivity  information.  Unfortunately  we  have  been  unable  to 
implement  the  direct  method  successfully  for  this  reaction  set,  presumably 
because  of  numerical  difficulties,  so  this  estimate  may  be  subject  to  some 
error.  The  Fourier  Amplitude  Sensitivity  Test, 11< 12  which  defines 
sensitivities  in  a  somewhat  different  fashion  by  averaging  over  the  uncertain¬ 
ties  of  other  parameters  in  the  model,  would  require  ~J0, 000  kinetic  runs 
for  a  57  parameter  system  for  fourth-order  interferences  according  to  a 
formula  given  by  its  authors.89 
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Detailed  Results  of  the  Ho-09  Calculations 


The  results  described  in  this  section  are  for  calculations  of  the  same 


type  as  that  described  in  the  previous  section.  That  is  to  say,  the  same 


error  tolerances,  integration  procedures,  etc.,  were  used,  and  similar 


numerical  behavior  was  observed.  Any  exceptions  are  specifically  noted.  In 


these  calculations  the  regions  around  each  of  the  three  explosion  limits  are 


given  particular  attention.  The  effects  of  inert  gas  diluents,  mixture 


composition,  temperature  and  pressure  on  the  behavior  in  the  vicinity  of 


these  limits  were  also  investigated.  Mechanistic  details  were  probed 


chiefly  by  sensitivity  analysis. 


As  a  preview  to  this  section,  consider  Figure  5,  which  displays  the 


results  of  several  kinetic  calculations  in  the  vicinity  of  the  three  limits. 


The  qualitative  and  some  quantitative  features  observed  in  experimental 


studies16-23  are  also  evident  in  these  modeling  results,  as  indicated  by  the 


trends  in  the  calculated  kinetic  behavior  shown  in  Figure  5 .  The  reproduction 


of  these  features  attests  to  the  validity  of  the  reaction  mechanism  and  the 


calculations . 


The  first  explosion  limit. 


In  the  region  just  above  and  below  the  first  limit  pressure,  diffusion 


appears  to  be  an  important  part  of  the  mechanism.  In  this  regard  the 


limitations  of  our  model  are  evident.  To  assess  the  effects  of  diffusion. 


calculations  are  often  carried  out  with  and  without  the  wall  reactions  (58  -62) 


The  assumption  is  that  wall  reactions  and  diffusional  mechanisms  achieve  a 


similar  result:  destruction  of  radical  species  before  the  reaction  is 


. 


Hill 
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completed.  The  general  validity  of  this  assumption  is  doubtful  (e.g.,  no 
account  is  allowed  for  the  surface- to-volume  ratio  in  the  wall  reaction 
model).  Nevertheless,  wall  reactions  are  included  to  test  their  effects. 

First  consider  a  calculation  carried  out  at  750°K  at  an  initial  pressure  - 
of  5.2  torr,  consisting  only  of  the  reactants  H2  and  02  present  in  a 
stoichiometric  ratio.  A  large  concentration  of  radicals,  particularly  H, 
develops  rather  quickly  after  a  short  induction  period  of  1  sec,  imnediately 
followed  by  an  explosion.  Sensitivities  of  [h]  to  the  rate  constants  were 
probed  at  0.15  seconds.  Here  and  in  further  discussions  of  sensitivities 
use  is  made  of  the  notation  a(b),  where  a  refers  to  the  reaction  number  of 
Table  I  whose  rate  constant  is  considered,  and  b  is  the  numerical  value  of 
the  relevant  normalized  sensitivity  coefficient.  The  sensitivities  are 
generally  rank-ordered  according  to  magnitude.  The  results  for  [h] 
sensitivities  are 


11(17)  >  20(1.4)  >  35(1.0)  >  12( -0 .7 )  >  26(0.1) 


Other  important  reactions  include  2,  3>  ^  and  24.  Reaction  35  1®  a  chain- 
initiator,  while  the  other  reactions  involve  chain-branching,  chain- 
propagation  and  chain-termination.  The  large  sensitivicy  of  Reaction  11 


(H  +  0£ 


OH  +  0)  indicates  the  importance  of  chain-branching  in  this  run. 


If  the  sensitivity  to  reaction  12  (H  +  02  +  M  — >  H02  +  M)  were  larger,  a 
slower,  less  explosive  reaction  would  be  indicated.  As  we  shall  see  later, 
the  branching  ratio  R  =  k.j/kioM,  obtained  from  a  simple  stability  theory, 

D 

is  an  Important  leading  indicator  of  whether  or  not  explosion  will  occur.21'70 

In  this  case  R  ~  1,  indicating  an  explosion.  When  the  calculations  are 
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repeated  Including  the  wall  reactions  (58-62),  the  overall  rate  Is  about  a 
factor  of  three  slower,  as  the  reactive  radical  species  are  destroyed. 

While  we  have  chosen  literature  values  of  the  wall  reactions  corresponding 
to  typical  situations,  it  is  true  that  certain  types  of  walls  can  promote 
radical  formation  rather  than  inhibiting  it.  Sensitivities  to  the  wall 
reactions  involving  U  and  0  (reactions  60  and  6l)  were  quite  large.  Other 
than  these  effects,  the  wall  reactions  change  the  major  features  of  the 
kinetic  behavior  very  little:  a  chain  mechanism  involving  virtually  the 
same  set  of  reactions  is  still  suggested  by  the  sensitivity  results.  A 
third  calculation  was  carried  out  with  the  same  initial  conditions  except 
that  700  torr  of  argon  was  present  as  a  diluent.  In  this  case  the  reaction 
rate  increases  quite  substantially.  This  effect  is  consistent  with  the 
experimental  result17*21  that  increased  inert  gas  pressure  causes  a 
lowering  of  the  first  explosion  limit  pressure  (initial  pressure  of  the  H2 
and  02.  reactants) .  This  pressure  lowering  is  generally  explained17  by  the 
argument  that  the  presence  of  inert  gas  hinders  diffusion  of  radicals  to  the 
wall  and  hence  promotes  rapid  reaction.  It  is  especially  interesting  that 
the  present  kinetic  calculations  indicate  an  increased  rate,  since  the  presence 
of  argon  might  be  expected  to  inhibit  a  branched-chain  explosion  due  to  the 
increase  in  [M]  which  decreases  R  . 
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The  effect  of  reactant  pressures  was  also  investigated.  At  750°K  with 
H2/02  present  in  a  stoichiometric  ratio  and  no  diluent  we  carried  out 
calculations  at  reactant  pressures  from  52  torr  down  to  £  0.1  torr.  At  the 
lowest  pressures  the  reaction  is  almost  three  orders  of  magnitude  slower 
than  for  the  highest  pressures.  Over  the  entire  pressure  range  essentially 
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the  same  set  of  reactions  were  important.  This  indicates  that  a  chain 
mechanism  persists,  yet  if  diffusion  were  included  in  the  model,  it  is  likely 
that  the  chain  would  never  propagate  to  any  appreciable  extent.  With  regard 
to  initial  pressures,  we  obtain  the  result  that  virtually  all  species 
concentrations  are  about  an  order  of  magnitude  more  sensitive  to  the  initial 
concentration  of  02  than  to  that  of  H2  In  this  regime.  This  result  likely 
reflects  (again!)  the  importance  of  reaction  11,  which  directly  involves  02 
but  not  H2. 

The  effect  of  temperature  was  also  studied.  As  indicated  by  the 
experiments, 17,21  the  first  limit  pressure  has  a  mild  temperature  dependence. 
The  calculations  indicated  the  first  limit  pressure  to  be  about  1  torr  or  so 
in  the  750  "  850°K  regime.  The  overall  reaction  rate  did  not  change  very 
much,  and  the  chain  mechanism  persists.  The  effects  of  including  wall 
reactions  in  the  mechanism  as  well  as  the  effects  of  the  presence  of  an  inert 
gas  were  similar  to  those  described  earlier.  On  the  other  hand,  at  slightly 
lower  temperatures,  particularly  where  the  first  and  second  limits  coalesce, 
some  dramatic  changes  are  observed.  For  example,  at  730°K  and  an  initial 
reactant  pressure  of  6  torr,  a  fast,  explosive  reaction  was  observed.  At 
650°K  for  the  6  torr  initial  reactant  pressure,  a  somewhat  slower,  yet  still 
explosive,  reaction  is  indicated.  However,  at  J20°K  and  an  initial  reactant 
pressure  of  14  torr,  a  very  slow  nonexplosive  reaction  results,  dominated  by 
the  presence  of  H02  and  H202,  rather  than  H,  0  and  OH. 


B.  The  second  explosion  limit. 


There  have  been  many  experimental  studies  o£  the  second  explosion  limit 
and  the  slow  reaction  between  the  second  and  third  limits.  Part  of  the 
reason  for  this  interest  is  that  chemical  kinetics  rather  than  fluid  dynamics  ~ 
is  expected  to  dominate  the  process.  Still,  the  quantitative  experimental 
results  obtained  for  the  limit  pressures  are  dependent  upon  the  vessel 
composition  and  the  method  of  preparation  of  the  reaction  vessel.  For  this 
reason  the  present  discussion  of  the  effects  of  temperature,  pressure  and 
inert  diluents  is  largely  qualitative. 

Several  calculations  were  performed  at  320  torr  for  a  stoichiometric 
ratio  of  reactants  from  temperatures  in  the  range  400-l600°K.  The  reaction 
is  negligibly  slow  at  room  temperature  for  this  pressure  regime.  Figure  3> 
which  shows  the  concentration  profiles  for  the  sample  calculation  of 
Section  III,  exhibits  the  features  typical  of  the  slow  reaction  characteristic 
of  the  temperature  range  500  -  800°K.  Within  a  very  short  time  [H02]  rises  to 
a  steady  concentration  of  about  IO10  molecules fen? .  The  log  [h20]  and  the 
log  [H202]  increase  steadily  and  slowly  as  nearly  linear  functions  of  log 
(time).  The  more  reactive  H,  0  and  OH  radicals  are  present  in  almost 
negligible  amounts.  Besides  the  concentrations,  the  sensitivity  results  for 
this  entire  pressure-temperature  regime  share  many  similar  characteristics; 
their  magnitudes  strongly  suggest  certain  mechanistic  details.  Reactions 
11  through  14,  which  have  already  been  recognized  as  important  in  an  earlier 
study,64  exhibit  the  largest  sensitivities.  In  general,  the  chain  is 
initiated  by  reaction  35>  propagated  by  reaction  11,  and  then  rapidly 
quenched  by  reaction  12.  Reaction  40  restores  the  chain  and  reactions  I3, 


14,  etc.,  promote  eventual  reaction.  Obviously,  H02  and  H202  are  important 
species  In  this  regime  and  rate  constants  for  reactions  Involving  them  should 
be  more  carefully  determined.  When  the  five  wall  reactions  are  included,  the 
sensitivities  to  reactions  58  and  59  are  very  large.  With  minor  exceptions,  - 
the  rate  constants  for  the  fifteen  reactions  of  the  Baldwin,  et  al.,2'7  set 
exhibit  the  largest  sensitivities.  This  result  is  reassuring  to  klnetlcists 
who  have  devoted  much  time  and  effort  to  understanding  the  most  critical 
elementary  steps  and  also  encouraging  since  these  results  demonstrate  the 
power  of  sensitivity  analysis  in  elucidating  mechanistic  details.  The 
differences  between  the  Baldwin  mechanism  and  that  indicated  by  the 
sensitivities  (see  Table  III)  are  very  slight.  The  chain- initiating 

reactions  H2  +  02  — *  20H  and  Ha  +  02  - >  H  +  H02  as  well  as  the  reaction 

0  +  H202  — *  H20  +  02  appear  to  be  more  important  than  the  reactions 
H  +  H02  ^  20H,  H  +  H02  *  H20  +  0,  and  0  +  H202  —■ >  OH  +  H02  of  the 

Baldwin  model.  Furthermore,  when  the  kinetic  calculation  was  rerun  with  only 
the  fifteen  reactions  of.  Table  III,  virtually  no  changes  in  kinetic  behavior 
were  observed.  The  concentration  profiles  for  this  run  are  shown  in  Figure  6. 

More  dramatic  effects  due  to  temperature  are  illustrated  in  Figures  Ja, 
b  and  c,  which  correspond  to  kinetic  runs  carried  out  at  825°K,  840°K  and 
875°K  respectively.  Figure  Ja  resembles  Figures  5  and  5  at  least  qualitatively. 
Noticeable  differences  are  a  several  hundredfold  increase  in  overall  reaction 
rate  and  substantial  increases  in  the  concentrations  of  the  more  reactive 
radical  species,  particularly  H.  The  branching  ratio  Rfi  =  kj^/kj^M  was 
about  0.5*  The  concentrations  and  sensitivity  profiles  are  similar  to  those 
for  the  T 00°K  case;  however,  the  normalized  sensitivities  to  reactions  11 
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and  12  are  nearly  equally  large  In  magnitude  but  opposite  in  sign.  At 
840°K  (see  Figure  7b)  a  slightly  delayed  explosion  occurs.  Before  0.40 
seconds  the  kinetic  behavior  is  similar  to  that  at  lower  temperatures,  yet 
around  0.42  seconds  extremely  rapid  concentration  changes  occur.  The  output 
from  the  kinetic  run  indicates  the  concentrations  of  B,  0,  OH  and  H20  to 
increase  by  about  two  orders  of  magnitude  within  10"3  sec.  At  short  times 
normalized  sensitivities  for  reaction  11  and  12  are  very  large.  At  the 
explosion  time  of  0.420  seconds,  all  normalized  sensitivities  are  huge,  many 
larger  than  a  million.  A  small  variation  in  virtually  any  rate  constant 
changes  the  explosion  time  (but  not  the  explosive  behavior)  for  this  process. 
The  largest  sensitivities  are  for  the  rate  constants  for  reactions  2,  3,  4, 

5,  8,  9,  10,  11,  12,  20,  23  and  26.  A  chain-branching  explosion  is  indicated 
at  this  temperature  and  pressure,  which  seems  to  be  close  to  the  point  where 
the  second  and  third  limit  pressures  coalesce.  At  875°K  (see  Figure  7c)  the 
explosion  occurs  almost  immediately .  The  branching  ratio  is  R  *»  0.6,  and 
[H02]  and  [H202]  never  get  very  large.  Sensitivities  are  not  nearly  as 
large  as  those  for  the  840°K  case,  but  at  10“3  seconds  the  sensitivities  of 
[H20]  to  the  following  rate  constants  are  largest: 

ll(-lo.4)  >  12(7.3)  >  55(-0.40)  >  5(-0.28)  ~  20(-0.28) 

>  l(-0 .22)  >  26(0.22)  >  24 (-0.16)  >  3(-0.l4)  >  2(0.11) 

The  mechanism  suggested  by  the  sensitivity  analysis  was  quite  similar  to  that 
suggested  by  screening  analysis.  This  was  often  the  case  for  many  of  the 
calculations  described  in  this  paper.  The  mechanism  suggested  by  the 
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sensitivities  for  this  regime,  however,  differs  considerably  from  that  of 
Baldwin,  et  al for  the  slow  reaction.  Since  little  attention  seems  to  have 
been  given  in  the  literature  to  the  mechanism  for  the  explosive  reaction 
regime,  we  have  developed  the  mechanism  given  in  Table  IV.  Note  especially  - 
the  importance  of  the  H  radical.  The  reaction  at  320  torr  for  temperatures 
greater  than  87  5°  K  exhibits  similar  explosive  characteristics.  Additionally, 
the  overall  reaction  rate  increases  and  normalized  sensitivities  generally 
decrease  slowly  as  the  temperature  increases. 

Several  calculations  were  performed  at  750°K  t0  probe  how  the  reaction 
proceeds  for  various  initial  pressures  of  the  reactants  in  the  neighborhood 
of  the  second  explosion  limit.  For  750°&  the  calculations  predict  an 
explosion  limit  pressure  of  80  ±  10  torr  for  a  stoichiometric  ratio  of  the 
reactants.  This  is  in  excellent  agreement  with  results  from  several 
experiments,  considering  the  effects  of  different  vessel  coatings  used.17'21'2^71 
The  mechanism  for  the  fast  reaction  below  this  limit  pressure  suggested  by  the 
analysis  of  sensitivities  is  exactly  that  given  in  Table  IV.  Above  the  limit, 
the  mechanism  for  the  slow  reaction  is  that  of  Table  III,  as  one  might 
reasonably  expect.  Perhaps  more  striking  is  that  fact  that  a  general  decrease 
in  overall  reaction  rate  is  observed  as  the  second  explosion  limit  pressure 
is  approached  from  the  high  pressure  side;  that  is,  the  reaction  is  much 
slower  at  100  torr  than  at  several  atmospheres.  This  is  in  semi>quantitative 
agreement  with  the  experiments  and  calculations  of  Lewis  and  Von  Elbe.21  A 
partial  explanation  is  the  dominance  of  reaction  12  suggested  by  its  large 
sensitivities  in  this  regime.  The  effect  of  the  concentration  of  M  in  the 
flux  term,  kj^MjCHlEOg]  ,  is  significant.  Although  just  before  the  limit. 
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the  reaction  is  very  slow,  at  reactant  pressures  below  80  torr,  the  overall 
reaction  rate  increases  several  orders  of  magnitude.  Large  sensitivities  to 
the  H  +  02  - *  OH  +  0  step  are  observed  and  the  branching  ratio  R  0.5. 

D 

An  "avalanche"  of  reactive  radicals  promotes  rapid  reaction. 

The  total  pressure  seems  to  play  as  significant  a  role  as  the  initial 
partial  pressures  of  H2  and  02.  Addition  of  50  -  100  torr  of  argon  or  nitrogen 
inhibits  the  explosion,  lowering  the  limit  pressure  to  below  30  torr  at 
750°K  for  the  sum  of  the  reactant  pressures .  This,  too,  is  in  agreement  with 
an  experiment.71  The  increased  concentration  of  M  makes  the  contribution 
to  reaction  12  more  dominant.  This  serves  to  slow  down  the  reaction 
considerably  just  above  the  limit,  but  eventually,  as  higher  pressures  are 
reached,  the  presence  of  argon  and/or  nitrogen  speeds  up  the  reaction  mildly. 
Several  third  body  reactions,  especially  reactions  12  and  Uo,  are  important 
in  the  formation  of  H20  via  the  H02  - *  H202  — *  H20  pathway. 

Several  other  miscellaneous  effects  on  the  kinetic  behavior  in  this 
regime  are  also  considered.  The  calculations  indicate  that  the  explosion 
limit  pressure  increases  steadily  with  temperatures  from  about  It  torr  at 
730°K,  where  the  first  and  second  limits  coalesce,  to  320  torr  at  840°K, 
where  the  second  and  third  limits  meet.  With  regard  to  deviations  from 
stoichiometric  ratios  of  reactants,  a  mild  increase  in  the  second  limit 
pressure  is  observed  for  an  increase  in  the  amount  of  hydrogen  in  the  mixture. 
Analysis  of  initial  concentration  sensitivities  indicates  that  most  of  the 
species  are  slightly  more  sensitive  to  the  initial  concentration  of  H2  than 
to  that  of  02.  This  is  difficult  to  explain,  although  the  dominance  of 
reaction  12,  H  +  02  +  M  — *  H02  +  M  ,  in  this  regime  certainly  plays  a 


role.  Perhaps  02  Is  less  important,  since  it  is  present  in  a  large, 
essentially  constant  concentration,  whereas  H  is  critical,  because  of 
reaction  12,  and  requires  H2  for  formation.  The  effects  of  temperature  and 


stoichiometry  are  in  agreement  with  experimental  results  ,17> 21*72*73 


C .  The  third  explosion  limit. 

The  region  in  the  vicinity  of  the  third  explosion  limit  has  been 
particularly  difficult  to  study  in  experiments.17*21’22  As  we  discussed  in 
Section  II,  both  chain  and  thermal  mechanisms  seem  to  play  a  role.  He  have 
carried  out  two  types  of  calculations  to  probe  this  regime;  one  type  assumes 
isothermal  behavior,  the  other  assumes  adiabatic  behavior.  It  is  likely  that 
the  actual  physical  situation  is  somewhere  between  these  two  extremes.  Heat 
and  mass  transport  undoubtedly  play  a  role  as  well. 

In  the  previous  subsection,  we  observed  that  an  increase  in  pressure 
corresponds  to  an  increase  in  overall  reaction  rate  in  between  the  second  and 
third  limits.  The  isothermal  calculations  indicate  that  this  behavior 

persists  up  to  at  least  one  hundred  atmospheres,  and  still  no  chain  branching 
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explosion  is  indicated.  Moreover,  no  sharp  limit  seems  to  exist  to  define 
the  third  limit.  Although  the  reaction  proceeds  very  rapidly,  concentration 
profiles  are  never  spiked,  but  remain  smooth  and  well-behaved.  The 
sensitivities  indicate  that  reactions  12,  1J,  14  and  40  -  the  same  as  of 
those  at  much  lower  pressures  -  exhibit  large  sensitivities  for  most  species. 
The  concentrations  of  H02  and  H202  are  orders  of  magnitude  higher  than  those 
of  H,  0  and  OH. 


The  adiabatic  calculations  are  probably  more  useful  in  this  regime.  For 
these  calculations  a  differential  equation  for  the  temperature  is  developed. 
(See  the  Appendix  for  a  detailed  derivation.)  Numerical  difficulties  were 
encountered  in  solving  the  combined  set  of  coupled  equations  [Eqs.(3.l)  and 
(A.6)].  The  temperature  changes  very  rapidly,  which  in  turn  can  cause  even 
more  rapid  changes  in  the  rate  constants  according  to  Eq.(2.l).  Other 
numerical  modelers  have  had  trouble  with  this  type  of  calculation.74  To 
remedy  this  situation,  we  generally  solved  a  difference  equation  for  T 
similar  to  Eq.(A.6)  at  each  step  taken  by  the  differential  equation  solver 
in  treating  the  usual  set  of  Eqs.(3«l)>  This  method  does  not  always  reproduce 
the  adiabatic  flame  temperature  (~5800°K),  which  is  defined  as  the  highest 
temperature  reached  assuming  that  all  the  internal  energy  is  used  to  raise 
the  temperature  of  the  reaction  products.  In  the  few  cases  where  we  solved 
the  differential  equations  [Eqs.(3.l)  and  (A.6)]  without  numerical  difficulty, 
the  results  obtained  agreed  with  those  obtained  by  the  difference  equation 
method.  The  difference  equation  method  should  give  an  accurate  temperature 
and  concentration  profile  for  at  least  the  Initial  1500°K  rise  in 
temperature . 

The  adiabatic  calculations  that  have  been  carried  out  in  between  the 
second  and  third  explosion  limits  share  certain  characteristics  in  common. 
First  of  all,  the  adiabatic  calculations  all  indicate  about  an  order  of 
magnitude  increase  in  overall  reaction  rate  in  comparison  to  the  isothermal 
calculations  at  the  same  initial  conditions.  In  each  run  there  is  an 
induction  period  which  may  be  very  long  in  the  case  of  the  reaction  just 
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above  the  second  limit  pressure  (~104  sec)  or  very  short  at  high  pressures 


and  temperatures  (&0.5  sec).  After  this  induction  period,  in  which 
the  temperature  rises  typically  a  few  hundred  degrees,  there  is  a  period  of 
explosive  reaction  accompanied  by  a  sharp  temperature  rise  of  a  few  thousand 
degrees.  At  these  temperatures  the  nascent  H20  quickly  reacts  to  form  H,  0 
and  OH  radicals,  which  are  relatively  more  stable  thermodynamically. 

Both  sensitivity  analysis  and  screening  analysis  were  used  to  elucidate 
mechanistic  details.  We  investigated  sensitivities  to  [A^}  ,  which  is 
temperature- independent,  rather  than  {k^}  .  The  analysis  of  sensitivities 
and  fluxes  indicates  that,  in  the  induction  period,  the  usual  mechanism  for 
the  reaction  between  the  second  and  third  limits  given  in  Section  IV. B  is 
most  appropriate.  On  the  other  hand,  careful  analysis  suggests  a  different 
mechanism  for  the  second  stage  of  rapid  reaction  and  sharp  temperature  rise. 
Radical  concentrations  of  H  and  0  increase  considerably  in  this  regime 
indicating  the  importance  of  several  chain  reactions,  most  notably  10,  11,  20, 
22,  2k,  25,  26  and  31.  At  such  high  temperatures  reactions  with  large 
activation  energies  like  5>  18,  19,  56  and  39  also  exhibit  reasonably  large 
sensitivities.  It  is  clear  that  many  elementary  reaction  steps  contribute  in 
the  vicinity  of  the  third  explosion  limit,  as  the  usual  slow  mechanism  is 
supplanted  by  a  chain  mechanism  including  even  some  endothermic  reactions. 

On  the  basis  of  these  calculations  a  few  conclusions  can  be  reached 
regarding  the  thermal  regime  studied  here.  First  there  does  not  seem  to  be 
any  clear-cut,  sharp  third  limit  pressure-temperature  boundary,  as  for  the 
first  and  second  limits.  The  calculations  do  indicate  that  at  high  pressures 
and  high  temperatures,  there  are  short  induction  periods  and  very  large 
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reaction  rates.  At  lower  pressures  and  temperatures,  induction  periods  are 


long  and  initial  overall  reaction  rates  are  small.  These  results  are  in 
qualitative  agreement  with  the  experimental  result  described  in  Section  II. 
According  to  several  theories,16'21  the  length  of  the  induction  period  and 
the  initial  reaction  rate  -  which  can  be  determined  for  various  pressures  and 
temperatures  from  calculations  -  are  key  Indicators  of  whether  an  explosion 
will  or  will  not  occur.  The  results  of  this  study  lend  credence  to  the  idea 
that  at  low  temperatures  and  pressures  induction  periods  are  so  long  that 
the  heat  may  be  transported  to  the  walls  of  the  reaction  vessel  and 
dissipated  before  reaction  occurs.  At  higher  pressures  and  temperatures 
reaction  occurs  before  the  heat  can  be  dissipated. 

As  an  epilogue  to  this  section  we  present,  in  Table  V,  a  list  of  the 
reactions  from  the  reaction  set  of  Table  I  which  seem  to  play  a  minor,  if 
any,  role  in  the  kinetics  of  the  process  regardless  of  the  temperature- 
pressure  regime  considered.  This  information  was  obtained  from  a  careful 
analysis  of  sensitivities.  Reruns  of  a  few  calculations  indicated  them  to 
be  of  very  little  significance.  Brief  explanations  for  their  unimportant 
stature  are  also  suggested. 


Sensitivity  Analysis  of  an  Experiment  to  Determine  the 
Rate  Constant  for  H  +  0P  +  H  — >  H0o  +  M 

We  describe  here  the  analysis  o£  a  shock  tube  experiment  carried  out  at 
temperatures  in  the  range  9$0  -  1176°K  by  Slack.40  In  this  study  experimentally 
determined  induction  times  were  matched  with  numerically  predicted  induction 
times  to  obtain  a  rate  constant  for  the  important  reaction  12  in  Table  I.  A 
model  containing  sixteen  reaction  steps  was  used  for  the  numerical  predictions. 
Slack  investigated  how  variations  in  the  numerical  values  of  rate  constants  in 
the  model  affected  the  rate  constant  for  H  +  02  +  M  — *  H02  +  M  (brute  force 
approach) .  Our  motivation  for  studying  this  experiment  was  to  demonstrate 
the  power  and  utility  of  derived  sensitivities  for  a  particularly  suitable 
experimental/theoretical  study  relevant  to  the  H2-02  system.  The  derived 
sensitivities,  which  have  already  been  considered  in  an  earlier  paper, 14 
were  compared  with  Slack's  results. 

To  simulate  one  of  the  shock  tube  experiments,  we  have  carried  out  two 
Isothermal  calculations  at  993°K.  The  initial  pressure  of  the  reactants, 
which  are  present  in  a  stoichiometric  ratio,  is  64  torr.  Nitrogen  (64  torr) 
was  used  as  a  diluent.  Three  calculations  were  carried  out:  the  first  used 
the  usual  57  reactions  of  Table  I,  the  second  included  only  those  sixteen 
reactions  in  Slack's  model  with  his  rate  constant  parameters,  and  the  third 
included  those  reactions  but  used  the  rate  constant  parameters  of  Table  I 
for  them.  No  appreciable  differences  are  detectable  for  the  three  runs. 

The  calculated  concentration  profiles  for  the  first  run  are  shown  in 
Figure  6.  The  calculated  induction  time  is  5  * 10"4  sec,  in  excellent 
agreement  with  the  experimental  result,  which,  according  to  the  oscilloscope 


trace  given  in  reference  40,  is  at  most  10#  higher.  Furthermore  the  net 
conversion  rates  (i.e.,  flux  contributions)  are  in  semi-quantitative 
agreement  with  those  given  by  Slack. 


The  derived  sensitivities  that  are  suitable  for  analyzing  this  experiment 
•re  of  the  form  (?jj“)  [or  ]  •  The  notation1*  I  refer, 

to  all  the  rate  constants  except  k^  (or  [H20])  and  k^  .  k^  is  determined 
in  this  experiment;  thus  it  is  indeed  appropriate  to  consider  it  an  output 
quantity.  [h20]  is  properly  considered  as  an  input  parameter,  since  its 


emission  spectrum  is  monitored  in  the  experiment.  The  relevant  dimensionless 
derived  sensitivity  coefficients  for  1 

U*"V&[h2o]  L  j 

have  been  calculated  and  rank-ordered.  The  results  at  4  X  10"4  are: 


11(26)  >  26(7-5)  >  20(4. 5)  >  [H20](2.4)  >  35(2.4)  >  2(-8.1XlO“2) 
>  24(-4.oxlO"2)  >  34(2.4  xlo"2)  >  4(1.5  xlo-2) 

£  3(1.5  Xlo;2)  >  43(-1.2x  lo"2)  >  i6(-3.6  x  lo"3) 


The  large  sensitivity  of  k ^  to  ku  is  expected.  The  experiment  is  carried 
out  in  high  temperature  regime,  where  reaction  11  dominates.  The  brute  force 
sensitivities  suggest  the  rate  constants  for  the  following  reactions  to  be 
important  in  determining  k12  :  4,  11,  13>  14,  20,  26  and  40.  Not  all  of  the 

results  of  reference  40  agree  with  ours.  Reactions  13,  14  and  40  -  important 
in  the  slow  regime  between  the  second  and  third  explosion  limits  -  are  not 
nearly  so  critical  for  the  rapid  reaction  characteristic  of  this  pressure- 

1 

temperature  regime.  On  the  other  hand,  the  chain- initiating  reaction,  35;  Is 


more  important  than  the  brute  force  results  indicate.  A  possible  reason 
for  the  discrepancies  observed  is  that  brute  force  sensitivity  analysis  is 
somewhat  difficult  to  Implement  systematically  even  for  small  models. 

Clearly,  to  determine  by  matching  experimentally  determined  and  numerically 
predicted  induction  times,  it  is  imperative  to  have  the  most  accurate  rate 
constant  parameters  available  for  reactions  11,  20,  26  and  35>  since  even 
slight  errors  can  introduce  large  uncertainties  in  because  of  the  high 


correlations 


VI.  Study  of  a  Related  Reaction:  The  Decomposition  of  Hydrogen  Peroxide 

To  demonstrate  the  versatility  o£  the  reaction  set  presented  in  Table  I, 
we  present  results  of  a  few  calculations  performed  to  study  the  decomposition 
of  hydrogen  peroxide,  a  process  which  exhibits  many  similarities  to  the  slow 
reaction  of  hydrogen  and  oxygen  between  the  second  and  third  limits.  The 
net  reaction  is 

H202(g)  - »  H20(g)  +  £02(g) 


Besides  this  process,  other  less-studied  reactions  (e.g.,  2H02  — *  H202  +  02) 
involved  in  the  H2-02  system  could  benefit  from  the  type  of  analysis  presented 
in  this  paper. 

Two  isothermal  calculations  were  carried  out  at  a  temperature  of  750°K. 
Partial  pressures  of  H202  and  argon  were  213  torr  and  360  torr,  respectively. 
For  the  first  calculation  the  usual  37  reaction  set  from  Table  I  was  utilized. 
The  resulting  concentration  profiles  are  given  in  Figure  9.  A  slow,  steady 
reaction  possessing  many  features  of  the  regime  between  the  second  and  third 
limits  for  hydrogen  combustion  is  indicated.  As  in  the  usual  slow  reaction, 
[H02]  rapidly  reaches  a  steady  concentration  and  remains  relatively  constant 
for  several  seconds.  A  previous  study75  assumed  [H02]  ~  0.0l[H202];  this 
apparently  is  only  true  after  several  hundred  seconds.  Virtually  no  chain¬ 
branching  is  indicated:  R_  £  0.1,  and  H  and  0  are  present  in  negligible 
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amounts.  Normalized  sensitivities  were  computed  for  several  times  and  all 
species  concentrations  for  this  reaction.  The  rank-ordered  normalized 
sensitivities  at  twenty  seconds  for  [OH]  are: 


pu- 


hO(l.O)  £  28(-1.0)  >  13(-1.8  X 10"3)  >  l6(-2.3  X 10-4) 


Most  of  the  other  rate  constants  exhibit  sensitivities  that  are  much  smaller. 
When  a  second  calculation  was  carried  out  including  only  the  four  reactions 
shown  above  as  the  model,  the  overall  rate  of  reaction  was  within  10$  of  the 
previously  mentioned  calculations,  although  some  intermediate  species 
concentrations  are  not  as  accurately  reproduced.  Several  other  reactions  are 
of  some  importance,  most  notably,  6,  7>  17,  18,  26  and  37*  Note  that  the 
mechanism  suggested  by  analysis  of  the  sensitivities  corresponds  roughly  to 
that  of  the  usual  slow  reaction  between  H2  and  02. 

Overall  the  agreement  with  relevant  experimental  results  and  earlier 
theoretical  models  is  reasonably  good.  Several  experiments  have  already 
established  the  nonchain  character  of  this  reaction.7®  Meyer,  et  al.,77 
have  studied  this  process  at  somewhat  higher  temperatures,  and  it  is 
interesting  that  their  suggested  mechanism  consists  of  the  reactions 

H202  +  M  - ►  20H  +  M 


OH  +  H202 


H02  +  H20 


2H02 


h2o2  +  02 


and  possibly 


OH  +  H02  - *•  H20  +  02  . 


This  is  exactly  the  mechanism  suggested  by  the  sensitivities.  Meyer,  et  al.. 
applied  the  steady  state  approximation76  to  [OH]  and  [H02]  in  the  first 


three  reactions  above  to  obtain  the  result 


=  -x-[h][h2o£] 


(6.1) 


where  X  -  2k4Q .  We  find  this  to  be  an  adequate  representation  on  the  basis 
of  the  results  presented  here,  with  the  exception  that  X  is  approximately 
I.5IC40  rather  than  2k40. 
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VII.  Discussion  and  Conclusions 


To  Che  besC  of  our  knowledge,  this  work  represents  Che  most  comprehensive 


modeling  study  available  for  Che  H2-02  system,  which  has  been  carefully 


studied  experimentally  and  theoretically  for  many  years.  In  many  cases  the 


modeling  results  could  be  anticipated  by  experimental  results  and  scientific 


intuition,  yet  these  calculations  still  serve  an  important  function  by 


amplifying  and  clarifying  certain  important  features  of  this  system. 


An  important  contribution  is  Table  I,  which  lists  a  useful  reaction  set 


and  appropriate  rate  constant  parameters.  Although  an  excellent  compilation 


has  already  been  provided  by  Baulch,  et  al .,29  for  this  system,  the  present 


evaluations  are  useful  since  many  of  these  rate  constant  parameters  have  been 


revised  and  updated  since  1971  and  have  not  yet  been  assembled  into  a  single 


source.79  Also,  the  fact  that  modeling  results  for  this  reaction  set  agree 


so  well  with  experiment,  attests  to  its  validity. 


A  second  important  facet  of  this  work  is  the  extensive  use  of  sensitivity 


analysis.  The  results  of  this  work  demonstrate  the  power  and  utility  of  a 


proper  analysis  of  sensitivity  information.  Sensitivity  analysis  was  chiefly 


applied  in  this  study  to  probe  mechanistic  behavior.  Although  it  is  true  that 


most  of  the  mechanistic  details  could  be  inferred  from  experimental  evidence 


by  careful  analysis,  determination  of  rate  constant  sensitivities  by  a  simple 


calculation  can  reduce  the  amount  of  work  considerably.  It  is  difficult  to 


determine  which  rate  constants  are  important  by  studying  of  a  list  of 


possible  reactions.  (The  mere  magnitude  of  the  rate  constants  is  often  a 


misleading  indicator.)  Sensitivity  aids  in  this  determination,  so  that  more 


complicated  analysis  may  not  be  warranted.  The  mechanisms  suggested  by  the 


'•.y.y-.-v'-.-v/ 
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sensitivities  given  in  Tables  III  and  IV  should  be  particularly  useful  for 
additional  model  studies.  Moreover,  Section  VI  demonstrates  the  use  of 
derived  sensitivities  in  interpreting  and  analyzing  the  results  of  kinetics 
experiments  for  determining  rate  constants.  By  establishing  the  correlations 
of  the  desired  rate  constant  to  other  rate  constants  in  the  system,  one  may 
obtain  more  realistic  error  bounds  for  the  experimental  determination. 
Kinetics  experiments  could  conceivably  be  designed  to  minimize  correlations 
with  the  desired  rate  constants. 

The  methodology  and  results  presented  in  this  paper  should  be  useful  for 
additional  studies  of  this  model  as  well  as  for  other  combustion  systems. 

The  methodology  can  certainly  be  extended  to  any  model  one  wishes  to  study. 
An  algorithm  for  numerical  modeling  and  sensitivity  analysis  has  been 
constructed  and  extensively  tested  and  is  easy  to  implement  for  Isothermal 
mass  action  kinetic  calculations.14* 80  Besides  the  methodology,  the  H2-02 
kinetic  results  from  the  present  or  similar  studies  should  be  useful  to  the 
experimental  kineticist  in  determining  which  elementary  steps  to  concentrate 
his  efforts  on.  We  found  that  many  reactions  Involving  H02,  particularly 
reaction  12,  are  important  in  the  regime  between  the  second  and  third 
explosion  limits.  These  reactions  have  been  extensively  studied,81*82  yet 
clearly  some  additional  work  is  called  for.  The  uncertainty  estimates 
given  in  Table  I  could  be  combined  with  the  sensitivity  results  to  assign 
each  reaction  a  sensitivity-uncertainty  index  as  a  guideline  for  useful 
experiments.83  The  H2-02  kinetic  results  may  be  useful  in  combustion  studies 
of  other  systems.  Westbrook,  et  al.,84  have  found  that  the  branching  ratio 
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of  H  +  02  — *  OH  +  0  and  H  +  02  +  M  — > 


H02  +  M  is  critical  in 


determining  the  length  of  the  induction  period  for  a  numerical  model  of 
methane  combustion.  The  important  elementary  steps  in  the  H2-02  system 
may  very  well  be  the  ones  that  are  critical  in  more  complicated  systems. 

We  conclude  this  work  with  a  few  suggestions  for  future  work  suggested 
by  our  results  on  this  and  similar  combustion  systems.  It  certainly 
would  be  useful  to  study  more  carefully  the  nature  of  explosions  character¬ 
istic  of  many  combustion  systems.  Obviously  explosions  result  from  some 
sort  of  instability.  For  the  H2-02  system  this  instability  is  evident  in 
the  calculated  concentration  profiles  (Figure  j)  and  in  the  sensitivity 
behavior  in  the  neighborhood  of  the  explosion  limits.  Thus  catastrophe 
theory85  or  stability  theory, 06-88  which  has  already  been  applied  to  the 
H2-02  system, 80  might  be  useful  approaches  to  this  problem,  possibly  in 
conjunction  with  sensitivity  analysis.  The  results  of  the  present  work 
could  be  used  to  provide  the  most  salient  features  of  the  explosive  system. 

Aside  from  general  studies  of  explosive  phenomena,  the  results  of  this 
work  suggest  that  additional  studies  of  the  first  and  third  limits  of  the 
H2-02  system  would  be  useful.  Mass  diffusion  should  properly  be  included 
into  the  model  for  studies  of  the  first  limit,  and  heat  transport  is  clearly 
important  for  studying  the  third  limit.  A  sensitivity  analysis  of  the 
resulting  partial  differential  equations  involved  in  such  studies  is 
somewhat  involved,  but  work  is  progressing  in  this  area  and  should  prove  to 
be  important.90  The  results  of  the  present  study  -  concentrations, 
sensitivity  profiles,  and  the  reaction  sets  of  Tables  111  and  IV  -  should  be 
useful  guides  for  treating  these  two  limits  as  well  as  other  pressure- 
temperature  regimes. 


Derivation  of  the  Expression  for 


dT 


in  Adiabatic  Modeling  Calculations 


UL  - - - 

The  £irst  lav  of  thermodynamics  (including  the  effects  of  chemical 
reactions)  is  given  by 


dU  =  dq  +  dw  +  £  (#)  dst 

1=1  1 


(A.l) 


where  U  denotes  the  internal  energy  of  the  system,  dq  refers  to  the 
differential  change  in  the  heat  added  to  the  system,  dw  is  the  differential 

,  and  s^  refers  to  the  number  of  moles 
of  species  i.  Since  adiabatic  conditions  are  assumed  and  since  the  reaction 
vessel  is  of  constant  volume,  dq  =  dw  =  0.  Moreover, 
dU  =  dH  -  d(PV)  =  dH  -  V  dp ,  so  that  Eq.(A.l)  becomes 


work  done  on  the  system 


(-Jwdv) 


40 I©4-!  -  vI(*r)d8 1 


(A.2) 


Dividing  through  by  V  and  noting  that  dU  -  ^  si^v  ®  we  obtain 

I  dI  -  I  (&0  dyl  '  1  (§-)  dyl  <*-5) 


where  y.  =  s  /V  is  the  usual  i-th  species  concentration  and  C  is  the 

Vi 

i-th  specific  heat  capacity  at  constant  volume.  Since  the  ideal  gas  law  is 
assumed  to  be  valid  for  All  species,  we  obtain 


L  Lhi  -  RTJ  d^i 


dT 


(A.4] 


where  is  the  specific  enthalpy  for  the  i-th  species.  Now  we 

formally  differentiate  with  respect  to  time  to  obtain 


dT 

dt 


I  °v.  *1 


(A.5: 


Since  h.  and  C  are  assumed  to  be  linear  functions  of  temperature 
i  v . 


Note  that  Eq.(A.6)  is  coupled  to  the  usual  set  of  equations  for  the  species 
concentrations,  Eqs.(^.l). 
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2  H02  — *  H202  +  02  1.8  (13  )  500  1.0 

H2  +  H02  - *  H202  +  H  6.0(11)  9,300  0.6 

ho2  +  h2o  — *  H202  +  oh  1.8  (13)  15,100  0.4 
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Table  I.  (continued) 


Table  1.  (continued) 
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Table  £.  (continued) 


Temperature  range:  300  -  2500°K.  Units  are  in  terms  of  moles,  eta3,  sec  and 
degrees  K.  Numbers  in  parentheses  denote  powers  of  ten.  Missing  entries 
are  all  zero. 

M  denotes  all  other  species  present  as  third  bodies.  The  numerical  values 
for  the  third-order  reactions  here  are  for  M  =  Argon  or  Nitrogen.  We  took 
numerical  values  for  these  rate  constants  for  M  =  H,  H2,  etc.  to  be  1.6  X 
those  shown  for  the  case  where  M  =  Argon.  Reference  29  was  used  as  a 
guide  for  assigning  these  third-body  efficiencies. 

Only  rough  estimates  of  the  rate  constant  uncertainties  are  given. 

In  several  cases  the  authors  made  slight  modifications  to  the  numerical 
values  of  the  rate  constant  parameters  given  in  the  references.  Often  a 
compromise  was  reached.  This  explains  the  reason  for  two  or  three 
citations  for  a  single' elementary  step. 

The  reaction  numbers  marked  with  an  asterisk  (*)  are  those  of  the 
Baldwin,  et  al.,  reaction  set  of  reference  27. 

The  wall  reaction  rate  constants  (58-62)  are  for  ken  spherical  vessels 
coated  with  boric  acid.  No  estimates  of  uncertainty  are  given,  as  the 
values  of  the  parameters  depend  critically  on  many  factors  (see  text). 


Table  11.  Coefficients  for  Computing  Specific  Heat  Capacities, 

[C  },  and  Specific  Enthalpies,  {h.},  for  the 
Vi  1 

H2-02  System 


Species 

c« 

V1 

c(i) 

hio) 

hW 

hi 

H 

2.981 

0 

52,580 

0.6778 

ho2 

8.566 

8.099  (-4) 

5,050 

-0.6415 

h2 

4.791 

6.361  (-4) 

0 

0 

02 

5.726 

5.512  (-4) 

0 

0 

h2o 

6.495 

1.586  (-5) 

-583,200 

-0.7126 

0 

2.970 

3.452  (-5) 

34,520 

-0.0100 

OH 

5-071 

4.965  (-4) 

9,650 

-0.4579 

h2o2 

8.465 

4.198  (-3) 

-33,500 

-0.7768 

9.706 

7.464  (-4) 

34,400 

0.0291 

Table  IV.  Mechanism  for  Che  Second  Explosion  Limit  Indicated  by  Sensitivity 


Analysis 


Reaction 


Number® 

Reaction 

Reason  for  Importance 

1 

2H  +  M  — >  H2  +  M 

H-containing 

2 

H  +  H02  — 9  H2  +  02 

H- containing 

3 

H  +  H02  - >  H20  +  0 

H-containing 

k 

H  +  H02  — 9  2  OH 

H-containing 

5 

H  +  H20  - 7  H2  +  OH 

H-containing 

8 

H  +  0  +  M  — 9  OH  +  M 

H-containing 

10 

H  +  OH  +  M  — 9  H20  +  M 

H-containing 

n 

H  +  02  — 9  OH  +  0 

chain-branching 

12 

H  +  02  +  M  — 9  H02  +  M 

chain  terminating 

20 

H2  +  0  - ►  H  +  OH 

chain-branching 

2k 

OH  +  0  - 9  H  +  02 

Reverse  of  Reaction  11 

26 

H2  +  OH  - *  H20  +  H 

H-containing 

35 

H2  +  02  - 9  2  OH 

chain- init iat ing 

a 


Reaction  numbers  correspond  to  those  given  in  Table  I. 


Table  V.  Insignificant  Reactions  for  the  H2-02  System  Indicated  by 


Sensitivity  Analysis 


Reaction 

Number® 

Reaction 

y+i 

Reason  for^mportance 

17 

OH  +  HOj,  - 9  VLs0z  +  0 

other  products  more  likely 

29 

2  OH  — 9  H  +  H02 

radical-radical 

30 

2  OH  - 9  H2  +  02 

radical-radical 

32 

20H  +  M  - 9  H202  +  M 

radical-radical 

33 

OE  +  M  - r  0  +  H  +  M 

endothermic 

38 

H2  +  M  - 9  2H  +  M 

endothermic 

42 

02  +  M  - >  2  0  +  H 

endothermic 

44 

02  +  OH  - 9  H02  +  0 

endothermic 

46 

H20  +  OH  - >  H2  +  H02 

endothermic 

48 

h2o  +  02  9  H202  +  0 

endothermic 

52 

2  02  - ►  0  +  O3 

endothermic 

54 

OH  +  02  - 9  H  +  O3 

endothermic 

a 


Reaction  numbers  correspond  to  those  given  in  Table  I 


Figure  1.  Illustration  of  the  variation  of  the  reaction  rate  as  a  function  of 
the  initial  pressure  of  H2  and  02.  Px,  P2,  and  P3  denote  the  first, 


second  and  third  explosion  limit  pressures,  respectively. 

T  =  750° K. 

Figure  2.  Explosion  limits  of  a  stoichiometric  H2-02  mixture  in  a  spherical 
KCl-coated  vessel  of  7.4  cm  diameter.  In  this  and  in  subsequent 
figures  the  explosion  limit  pressure  is  the  sum  of  the  initial 
partial  pressures  of  H2  and  02  which  cause  explosion  at  a 
particular  temperature.  No  diluent  is  present.  The  dotted  lines 
indicate  that  the  boundaries  are  uncertain.  (After  reference  21, 

P.  9.) 

Figure  3.  Concentration  profiles  for  the  calculation  described  in  Section  III. 

T  =  700°K,  P^nit  =  320  torr,  stoichiometric  mixture  of  H2  and  02, 
no  diluent  present.  This  figure,  as  well  as  Figures  6,  7,  8  and  9 
were  produced  directly  from  the  differential  equation  solver  mesh 
points.  Straight  lines  are  drawn  between  the  irregularly  spaced 
points . 

Figure  4.  Plot  of  a  typical  integrand  for  a  sensitivity  coefficient  for  the 
sample  calculation  described  in  Section  III.  This  plot  is  for 
the  coefficient  3[H]/dk20  evaluated  at  103  sec.  The  scale  is 
expanded  for  the  final  2f>  of  the  interval.  The  shaded  area  under 

i 

the  curve  illustrates  the  contribution  to  the  Integral. 


Figure  5.  Same  aa  Figure  2,  except  that  the  result8  of  a  number  of  isothermal 


kinetic  calculations  are  displayed  as  well.  The  notation  is  as 
follows:  S  denotes  slow  reactions,  which  according  to  the 
calculations  take  more  than  10  seconds  to  achieve  50J&  completion;  . 
M  refers  to  moderately  fast  reactions,  which  take  between  0.5  and 
10  seconds  to  reach  50$  completion;  and  F  refers  to  fast  reactions, 
which  require  fewer  than  0 .5  seconds  to  reach  50 ^  completion.  The 
notation  R  refers  to  situations  for  which  large  concentrations  of 
reactive  radicals  H,  0  and  OH  are  observed.  Such  reactions  are 
generally  explosive.  N  refers  to  smaller  concentrations  of  the 
reactive  radicals  and  larger  concentrations  of  the  species  H02 
and  H202.  Such  reactions  are  generally  slower  (S  or  M)  and  display 
non-explosive  behavior. 

Figure  6.  Same  as  Figure  3,  except  that  only  the  15  most  sensitive  reactions 
are  used  as  the  model.  03  does  not  appear  since  it  does  not  enter 
into  any  of  these  reactions. 

Figure  7*  (a)  Concentration  profiles  for  a  calculation  carried  out  at 

T  =  825°K,  =  520  torr.  Stoichiometric  mixture  of  Hg 

and  02,  no  diluent  present.  Non-explosive  behavior  shown. 

(b)  Same  as  (a),  except  that  T  =  840°K.  A  delayed  explosion  is 
indicated. 

(c)  Same  as  (a),  except  that  T  =  87 5° K .  Explosive  behavior  is 
indicated . 


Figure  8.  Concentration  profiles  for  a  calculation  carried  out  at 

T  =  993°K,  p  *  64  torr.  Nitrogen  (64  torr)  is  present  as 
a  diluent.  This  calculation  is  for  the  experiment  described  in 
reference  40.  The  induction  time  is  obtained  from  the  rapid 
decrease  in  H2  pressure;  our  estimate  is  5xl0“4sec. 


Figure  9.  Concentration  profiles  for  the  decomposition  of  H202.  T  =  750R, 
Pinit(H  0  )  *  torr*  Argon  (560  torr)  is  present  as  a 


diluent.  The  concentrations  of  H,  0  and  O3  are  very  small  and  off 
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